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I discuss the impact of color-octet contributions and higher-order QCD corrections on 
the cross section for inelastic J/ip photoproduction. The theoretical predictions are 
compared with recent experimental data obtained at HERA. 

1. Introduction 

The production of heavy quarkonium states in high-energy colhsions provides an 
important tool to study the interplay between perturbative and non-perturbative 
QCD dynamics. While the creation of heavy quarks in a hard scattering process 
can be calculated in perturbative QCD0, the subsequent transition to a physical 
bound state introduces non-perturbative aspects. A rigorous framework for treating 
quarkonium production and decays has recently been dfiveloped.0 The factorization 
approach is based on the use of non-relativistic QCDB (NRQCD) to separate the 
short-distance parts from the long-distance matrix elements and explicitly takes 
into account the complete structure of the quarkonium Fock space. This formalism 
implies that so-called color-octet processes, in which the heavy-quark antiquark pair 
is produced at short distances in a color-octet state and subsequently evolves non- 
perturbatively into a physicBLquarkonium, should contribute to the cross section. 
It has recently been arguedau that quarkonium production in hadronic collisions 
at the Tevatron can be accounted for by including color-octet processes and by 
adjusting the unknown long-distance color-octet matrix elements to fit the data. 

In order to establish the phenomenological significance of the color-octet mecha- 
nism it is necessary to identify color-octet contributions in different production pro- 
cesses. Color-octet pxoduction of J/tp particles has also been studied in the context 
of e~^e~ annihilfitionEl, Z decayaJ, hadronic collisions at fixed-target experimcntsBu 
and B decaysE2l. Here, I review the impact of color-octet contributions and higher- 
order QCD corrections on the cross section for J/tp photoproduction. The produc- 
tion of J/tp particles in photon-proton collisions proceeds predominantly through 
photon- gluon fusion. Elastic/diffractive mechanismsllil can be eliminated by mea- 
suring the J/tJj energy spectrum, described by the scaling variable z = p ■ / p ■ k-y, 
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with p, k^_y being the momenta of the proton and J/tpi 7 particles, respectively. 
In the proton rest frame, z is the ratio of the J/ip to 7 energy, z = E^/Ej. For 
elastic/diffractive events z is close to one; a clean sample of inelastic events can be 
obtained in the range z < 0.9. 

According to the NRQCD factorization formalism , the inclusive cross section 
for J/ip photoproduction can be expressed as a sum of terms, each of which factors 
into a short-distance coefficient and a long-distance matrix element: 

dCT(7 -\-g^ J/i; + X)^J2 "^^(^ + g^cc + X) {O-'/^ [71]) (1) 

n 

Here, da denotes the short-distance cross section for producing an on-shell cc-pair 
in a color, spin and angular-momentum state labelled by n. The NRQCD matrix 
elements (O'^/'^ [n]) = (0|O'^/'^ [n]|0) give the probability for a cc-pair in the state 
n to form the J/ip particle. The relative importance oLthe various terms in (|l|) 
can be estimated by using NRQCD velocity scaling rules.Ej For v ^ {v being the 
average velocity of the charm quark in the J/ip rest frame) each of the NRQCD 
matrix elements scales with a definite power of v and the general expression (^) can 
be organized into an expansion in powers of f^. 

2. Color-singlet contribution 

At leading order in-U^ , eq. reduces to the standard factorization formula of the 
color-singlet modelliS. The short-distance cross section is given by the subprocess 

l + g^cc[h^Si]+g (2) 

shown in Fig.|l|a, with cc in a color-singlet state (denoted by 1), zero relative velocity, 
and spin/angular- momentum quantum numbers '^^^^ Lj — ^ Si. Up to corrections of 
0{v'^), the color-singlet NRQCD matrix element is related to the J/ip wave function 
at the origin through {O-^l'l' [h^Si]) « (9/27r)|(^(0)|2 and can be extracted from the 
measurement of the J/ip leptonic decay width or calculated within potential models. 
Relativistic corrections due to the motion of the charm quarks in the J/ip liDund 
state enhance the large- 2; region, but can be neglected in the inelastic domain.La The 
calculation of the higher-order perturbative QGILijorrections to the short-distance 
cross section (|^) has been performed recently.ll3't3 Generic diagrams which build 
up the cross section in next-to-leading order (NLO) are depicted in Fig.|l|. Besides 
the usual self-energy diagrams and vertex corrections for photons and gluons (b), 
one encounters box diagrams (c), the splitting of the final-state gluon into gluon 
and light quark-antiquark pairs, as well as diagrams renormalizing the initial state 
parton densities (e) . Inclusion of the NLO corrections reduces the scale dependence 
of the theoretical prediction and increases the cross section significantly, depending 
in detail on the 7p energy and the choice of parameters.Ej Details of the calculation 
and a comprehensive analysis of total cross sections and differential distributions 
for the energy range of the fixed-target experiments and for J/ip photoproduction 
at HERA can be found elsewhere.E3 
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Fig. 1. Generic diagrams for J/ip photoproduction via the color-singlet channel: (a) leading order 
contribution; (b) vertex corrections; (c) box diagrams; (d) splitting of the final state gluon into 
gluon or light quark-antiquark pairs; (e) diagrams renormalizing the initial-state parton densities. 



3. Color-octet contributions 

Color-octet conf igu rations are produced at leading order in as through the 2 — > 1 
parton processesOEllli 

7 + 5 ^ cciS/S-o] 

7 + 5 -> cc[8,3Po,2] (3) 

shown in Fig.||a. Due to kinematical constraints, the leading color-octet terms will 
only contribute to the upper endpoint of the J/ip energy spectrum, z « 1 and 
p± ~ 0, p± being the J/ip transverse momentum. Color-octet configurations which 
contribute to inelastic J/?i photoproduction z < 0.9 and p± > 1 GeV are produced 
through the subprocessesEZlE3 

1 + 9 cc[8, ^S'ol+g 
1 + 9 cc[8,'^Si]+g 

1 + 9 cc[8,3Po,i.2] +5 (4) 

as shown in FigJ^b. Light-quark initiated contributions are strongly suppressed at 
HERA energies and can safely be neglected. 
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(bi) (b.) 

Fig. 2. Generic diagrams for J/^p photoproduction via color-octet channels: (a) leading color-octet 
contributions; (b) color-octet contributions to inelastic J/ip production. 

The transition of the color-octet cc [8, ■^'^+^Lj] pair into a physical J/4' state 
through the emission of non-perturbative gluons is described by the long-distance 
matrix elements {O'-^^"^ They have to be obtained from lattice simula- 

tionscj or measured directly in some production process. According to the velocity 
scaling rules of NRQCD, the color-octet matrix elements associated with S'-wave 
quarkonia should be suppressed by a factor of compared to the leading color- 
singlet matrix element. ^Color-octet contributions to J/tp photoproduction can thus 
become important only if the corresponding short-distance cross sections are en- 
hanced as compared to the color-singlet nmcess. Color-octet matrix elements have 
been fitted to prompt J/V' data from CDFej and found to be O(10~^ GeV'^), consis- 
tent with the NRQCD velocity scaling rulesD^B Meanwhile, fit values for color-octet 
matrix elements have also been obtained from analyses of quarkonium production in 
hadrnnic collisions at fixed-target experirnentsO, J/ip photoproduction at the elastic 
peaUla and J/^J production in B decaysll2l. The results seem to indicate that the 
values for the color-octet matrix elements extracted from the Tevatron data at mod- 
erate p± are too large; they should however be considered with some caution since 
significant higher-twist corrections are expected to contribute in the small-pj^ region 
probed at fixed target experiments and in elastic J/ip photoproduction. Moreover, 
the comparison between the different analyses is rendered difficult by the fact that 
the color-octet matrix elements can in general only be extracted in certain linear 
combinations which depend on the reaction under consideration, see Sec. 4. 

4. J/ip photoproduction at HERA 

The production of J/ip particles in high energy ep collisions at HERA is domi- 

^In the case of i?-wave quarkonia, color-singlet and color-octet matrix elements contribute at the 
same order in v.tH Photoproduction of P-wa\j O ntc ttes is, however, suppressed compared with J/tp 
states, by two orders of magnitude at HERAJ23Ej 



Color- Singlet and Color- Octet Contributions to J/ip Photoproduction 5 



nated by photoproduction events where the electron is scattered by a small angle 
producing photons of almost zero virtuality. The measurements at HERA provide 
information on the dynamics of J/ip photoproduction in a wide kinematical re- 
gion, 30 GeV < y^s^p < 200 GeV, corresponding to initial photon energies in a 
fixed-target experiment of 450 GeV < < 20,000 GeV. Due to kinematical con- 
straints, the leading color-octet processes (|^) contribute only to the upper endpoint 
of the J/'0 energy spectrum, z « 1 and p± « 0. The color-singlet and color-octet 
predictions (||) have been compared to experimental datacJ obtained in the region 
z > 0.95 and p± < I GeV.EZi Since the factorization approach cannot be used to 
describe the exclusive elastic channel 7 -I- p — *■ J/ip + P, elastic contributions had 
been subtracted from the data sample. It was shown that the large cross section 
predicted by using color-octet matrix elements as extracted from the Tevatron fits 
appears to be in conflict with the experimental data. It is, however, difficult to 
put strong upper limits for the octet terms from a measurement of the total cross 
section in the region z 1 and p±^ « since the overall normalization of the the- 
oretical prediction depends strongly on the choice for the charm quark mass and 
the QCD coupling. Moreover, diffractive production mechanisms which cannot be 
calculated within perturbative QCD might contaminate the region z « 1 and make 
it difficult to extract precise information on the color-octet contributions. Finally, 
it has been argued that sizable higher-twist effects are expected to contribute in the 
region p± < 1 GeV, which cause the breakdown of the factorization formula (|l|) .c3 
It is therefore more appropriate to study J/ip photoproduction in the inelastic 
region z < 0.9 and p± > 1 GeV where no diffractive channels contribute and where 
the general factorization formula (0) and perturbative QCD calculations should be 
applicable. Adopting the NRQCD matrix elements as extracted from the fits to 
prompt J/ip data at the Tevatron one finds that color-octet and color-singlet con- 
tributions to the inelastic cross section are predicted to be of comparable size.0lli 
The short-distance factors of the [8, -^^o] and [8, •^Po,2] channels are strongly en- 
hanced as compared to the color-singlet term and partly compensate the 0(10^^) 
suppression of the corresponding non-perturbative matrix elements. In contrast, 
the contributions from the [8,^ Si] and [S,^Pi] states are suppressed by more than 
one order of magnitude. Since color-octet and color-singlet processes contribute at 
the same order in as, the large size of the [8, ^^o] and [8, ^^0,2] cross sections could 
not have been anticipated from naive power counting and demojistrates the crucial 
dynamical role played by the bound state quantum numbers .Ell As for the total 
inelastic cross section, the linear combination of the color-octet matrix elements 
(^QJ/i: [8, ^Sq]) and (O'^/'^ [8,^Po]) that is probed at HERA is almost identical to 
that extracted from the Tevatron fits at moderate p±, independent of y/s^.^ The 
Tevatron results can thus be used to make predictions for color-octet contributions 
to the total inelastic J/ip photoproduction cross section without further ambigui- 
ties. However, taking into account the uncertainty due to the value of the charm 



tAt leading order in v^, the P-wave matrix elements are related by heavy-quark spin symmetry, 
(CU/V- [8, 'ipj]) ^ (2 J + 1) (O-^/* [8, 3Po]>. 
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quark mass and the strong coupling, the significance of color-octet contributions 
cannot be deduced from the analysis of the absolute J/ip production rates. In 
fact, the experimental data can be accounted for by the color-singlet channel alone, 
once higher-order QCD corrections are included and the theoretical uncertainties 
due to variation of the charm quark mass and the strong coupling are taken into 
account, as demonstrated at the end of this section. The same statement holds 
true for the transverse momentum spectrum, since, at small and moderate p± , both 
color-singlet and color-octet contributions are almost identical in shape. At large 
transverse momenta, p± > lQ_GfiV, charm quark fragmentation dominates over the 
photon- gluon fusion process In contrast to what was found at the Tevatroni^l, 
gluon fragmentation into color-octet states is suppressed over the whole range of 
Pj^ in the inelastic region z < 0.9E3 

A distinctive signal for color<)ctet processes should, however, be visible in the 
J/ijj energy distribution dcr/dz.EII The linear combination of color-octet matrix el- 
ements that is probed by the J/iJj energy distribution does, however, depend on 
the value of z. Therefore, one cannot directly use the Tevatron fits but has to 
allow the individual color-octet matrix elements to vary in certain ranges, con- 
strained by the value extracted for the linear combination. It has in fact been 
argued that the color-octet matrix element {O'^^^AB, ^Po]) could be negative due to 
the subtraction of power ultraviolett divergences .Eil In contrast, the matrix element 
^Qj/ip jg^ ^5o]) is free of power divergences and its value is thus always positive. Ac- 
cordingly, I have allowed (O"'/'^ [8, ^Pq]) /ml to vary in the range [-0.01, 0.01] GeV^ 
and determined the value of the matrix element (O"^/'^ [8, ^^o]) from the linear 
combination extracted at the Tevatron. § The result is shown in Fig.||(a) where I 
have plotted (leading-order) color-singlet and color-octet contributions at a typical 
HERA energy of ^s^p — 100 GeV inthe restricted range p^>l GeV, compared 
to recent experimental data from HlEa and preliminary data from ZEUSO. The 
hatched error band indicates how much the color-octet cross section is altered if 
l^Qj/i' [8,3P(,])/m2 varies in the range [-0.01,0.01] GeV^ where the lower bound 
corresponds to {O-'/'^ %^Po\)/ml = -0.01 GeV^. Since the shape of the distri- 
bution is almost insensitive to higher-order QCD corrections or to the uncertainty 
induced by the choice for rric and as, the analysis of the J /ip energy spectrum dcr/dz 
should provide a clean test for the underlying production mechanism. From Fig.^ 
one can conclude that the shape predicted by the color-octet contributions is not 
supported by the experimental data. The discrepancy with the data can only be re- 
moved when reducing the relative weight of the color-octet contributions by at least 
a factor of five.EZi Let me emphasize that the rise of the cross section towards large 
z predicted by the color-octet mechanism is not sensitive to the small-pj^ region 
and thus not affected by the coUinear divergences which show up at the endpoint 
z = 1 and p±^ = 0. This is demonstrated in FigJ^(b) where I show color-singlet and 
color-octet contributions to the J /"ip energy distribution for pj^ > 5 GeV. It will be 

§Note that, given (£)■-'/''' [8, ^So]) < 0.1 GeV"^ as required by the velocity scaUng rules, a value 
(O^/'/' [8, ^Po])/m^ < - 0.01 GeV-* would bo in contradiction with the Tevatron fits. 
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Fig. 3. Color-singlet and color-octet contributions to the J/(/> energy distribution da/dz at the 
photon-proton centre of mass energy ^s^p = 100 (j oV i ntegrated in the range (a) > 1 GeV 
and (b) pj_ > 5 GeV compared to experimental dataE3c3- 

very interesting to compare these predictions with data to be expected in the future 
at HERA. Let me finally mention that the shape of the J/4> energy distribution 
couldJpe influenced by the emission of soft gluons from the intermediate color-octet 
state.tl While this effect, which cannot be predicted within the NRQCD factoriza- 
tion approach, might be significant at the elastic peak, it is by no means clear if and 
in which way it could affect the inelastic region z < 0.9 and p± > I GeV. In fact, 
if soft gluon emission were important, it should also result in a feed-down of the 
leading color-octet contributions (^) into the inelastic domain, thereby increasing 
the discrepancy between the color-octet cross section and the data in the large-z 
region. 

For the remainder of this section, I will demonstrate that the experimental 
results on differential distributions and total cross sections are well accounted for 
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by the color-singlet channel alone including higher-order QCD corrections. This can 
e.g. be inferred from Figj^ where I compare the NLO color-singlet predi cti on for 
the J/i/' transverse momentum distributionts with recent results from HlEj. Note 
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Fig. 4. LO and NLO color-singlet prediction for the J/tf> transverse momentum spectrum dcr/dp^ 
at the photon-proton centre ofp«ass energy y/s^yp = 100 GeV integrated in the range z < 0.9 
compared to experimental datacj. 

that the inclusion of higher-order QCD corrections is crucial to describe the shape 
of the pj_ distribution. However, a detailed analysis of the transverse momentum 
spectrum reveals that the fixed-order perturbative QCD calculation is not under 
proper control in the limit p± — > 0, Fig.^. No reliable prediction can be made in the 
small-p^ domain without resummation of large logarithmic corrections caused by 
multiple gluon emission. If the region p_L < 1 GeV is excluded from the analysis, the 
next-to-leading order color-singlet prediction accounts for the energy dependence of 
the cross section and for the overall normalization, Fig. ||. The sensitivity of the 
prediction to the small- a; behaviour of the gluon distribution is however not very 
distinctive, since the average momentum fraction of the partons < a; > is shifted to 
larger values when excluding the small-pj^ region. 

5. Conclusion 

I have discussed color-singlet and color-octet contributions to the production of 
J/^/i particles in photon-proton collisions, including higher-order QCD corrections 
to the color-singlet channel. A comparison with photoproduction data obtained at 
fixed-target experimentaiB and the ep collider HERA reveals that the J/^ energy 
spectrum and the slope of the transverse momentum distribution are adequately 
accounted for by the next-to-leading order color-singlet prediction in the inelastic 
region p± > 1 GeV and z < 0.9. Taking into account the uncertainty due to 
variation of the charm quark mass and the strong coupling, one can conclude that 
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Fig. 5. NLO color-singlet prediction for the total inelastic J/0 photp«roduction cross section as 
a function of the photon-proton energy for diSenfeint parametrizationsH of the parton distribution 
in the proton compared to experimental dataEjE3- 



the nomialization too appears to be under semi-quantitative control. Higher-twist 
effectso must be included to improve the quality of the theoretical analysis further. 
Distinctive signatures for color-octet processes should be visible in the shape of 
the J/'ijj energy distribution. However, these predictions appear at variance with 
recent experimental data obtained at HERA indicating that the values of the color- 
octet matrix elements {O'^^"^' [S,^ Sq]) and (O'^'^''^ [8, ■^Pq]) are considerably smaller 
than suggested by the fits to Tevatron data at moderate p±. Support is added to 
this result by recent analyses J/-0 production in hadronic collisions at fixed- 
target energieaJ and B decaysEiJ. Clearly, much more effort, both theoretical and 
experimental, is needed to establish the phenomenological significance of color-octet 
contributions to J/ip production and to proof the applicability of the NRQCD 
factorization approach to charmonium production in hadronic collisions at moderate 
transverse momentum. 
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